Comparing modern identification methods for wild bees: Metabarcoding and image-based morphological taxonomic assignment.

With the decline of bee populations worldwide, studies determining current wild bee distributions and diversity are increasingly important. Wild bee identification is often completed by experienced taxonomists or by genetic analysis. The current study was designed to compare two methods of identification including: (1) morphological identification by experienced taxonomists using images of field-collected wild bees and (2) genetic analysis of composite bee legs (multiple taxa) using metabarcoding. Bees were collected from conservation grasslands in eastern Iowa in summer 2019 and identified to the lowest taxonomic unit using both methods. Sanger sequencing of individual wild bee legs was used as a positive control for metabarcoding. Morphological identification of bees using images resulted in 36 unique taxa among 22 genera, and >80% of Bombus specimens were identified to species. Metabarcoding was limited to genus-level assignments among 18 genera but resolved some morphologically similar genera. Metabarcoding did not consistently detect all genera in the composite samples, including kleptoparasitic bees. Sanger sequencing showed similar presence or absence detection results as metabarcoding but provided species-level identifications for cryptic species (i.e., Lasioglossum). Genus-specific detections were more frequent with morphological identification than metabarcoding, but certain genera such as Ceratina and Halictus were identified equally well with metabarcoding and morphology. Genera with proportionately less tissue in a composite sample were less likely to be detected using metabarcoding. Image-based methods were limited by image quality and visible morphological features, while genetic methods were limited by databases, primers, and amplification at target loci. This study shows how an image-based identification method compares with genetic techniques, and how in combination, the methods provide valuable genus- and species-level information for wild bees while preserving tissue for other analyses. These methods could be improved and transferred to a field setting to advance our understanding of wild bee distributions and to expedite conservation research.

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population.

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

Influencia potencial de la divergencia adaptativa a nivel genoma sobre el resultado de la reubicación para conservación en una población aislada de urogallo mayor Resumen Las reubicaciones para conservación son una herramienta importante que se usa con frecuencia para aumentar las poblaciones en declinación o reestablecer las poblaciones erradicadas. Una de las metas de este aumento es incrementar la diversidad genética y reducir el riesgo de depresión endogámica (es decir, rescate genético). Sin embargo, la introducción de individuos de una población con divergencia significativa puede perturbar los rasgos coadaptados y reducir la aptitud local (es decir, depresión exogámica). La información genética es cada vez más accesible para las especies silvestres y puede proporcionar conocimiento único con respecto a la presencia y retención de la variación genética introducida a partir del aumento como un indicador de eficiencia y las similitudes adaptativas como un indicador de la idoneidad de la población de origen y la receptora. Usamos dos conjuntos de datos genéticos para evaluar el aumento de las poblaciones aisladas del urogallo mayor (Centrocercus urophasianus) en la región noroeste de la distribución de la especie (Washington, EUA) y para evaluar de forma retrospectiva la divergencia adaptativa entre la población de origen y la receptora. Desarrollamos dos modelos estadísticos para los datos microsatelitales para así evaluar los resultados del aumento. Usamos un modelo para predecir la diversidad genética después del aumento y comparamos estas predicciones con observaciones del cambio genético. Usamos el segundo modelo para cuantificar el aumento de la reproducción observada atribuida a las reubicaciones (evidencia de la integración poblacional). También caracterizamos la divergencia adaptativa a nivel genoma entre la población de origen y la población receptora. La diversidad genética observada (HO = 0.65) fue mayor de lo que se predijo en la población receptora de no haber ocurrido el aumento (HO = 0.58) pero menor de lo que se predijo en nuestro modelo (HO = 0.75). El aumento de la variación genética compartida entre las dos poblaciones residentes geográficamente aisladas incrementó, lo cual es evidencia de un flujo génico periódico que antes se supuso casi no ocurría. Entre los genes adaptativos candidatos asociados a una FST elevada (143 genes) o a variables ambientales locales (97 y 157 genes para cada método de asociación entre el ambiente y el genotipo, respectivamente) encontramos grupos de genes con funciones relacionadas que pueden influir sobre la habilidad de cada reubicación para usar recursos locales y navegar ambientes desconocidos y su potencial reproductivo, todas posibles razones para la baja retención genética en el aumento.

Assessing arthropod diversity metrics derived from stream environmental DNA: spatiotemporal variation and paired comparisons with manual sampling.

Background: Benthic invertebrate (BI) surveys have been widely used to characterize freshwater environmental quality but can be challenging to implement at desired spatial scales and frequency. Environmental DNA (eDNA) allows an alternative BI survey approach, one that can potentially be implemented more rapidly and cheaply than traditional methods. Methods: We evaluated eDNA analogs of BI metrics in the Potomac River watershed of the eastern United States. We first compared arthropod diversity detected with primers targeting mitochondrial 16S (mt16S) and cytochrome c oxidase 1 (cox1 or COI) loci to that detected by manual surveys conducted in parallel. We then evaluated spatial and temporal variation in arthropod diversity metrics with repeated sampling in three focal parks. We also investigated technical factors such as filter type used to capture eDNA and PCR inhibition treatment. Results: Our results indicate that genus-level assessment of eDNA compositions is achievable at both loci with modest technical noise, although database gaps remain substantial at mt16S for regional taxa. While the specific taxa identified by eDNA did not strongly overlap with paired manual surveys, some metrics derived from eDNA compositions were rank-correlated with previously derived biological indices of environmental quality. Repeated sampling revealed statistical differences between high- and low-quality sites based on taxonomic diversity, functional diversity, and tolerance scores weighted by taxon proportions in transformed counts. We conclude that eDNA compositions are efficient and informative of stream condition. Further development and validation of scoring schemes analogous to commonly used biological indices should allow increased application of the approach to management needs.

A genetic warning system for a hierarchically structured wildlife monitoring framework.

Genetic variation is a well-known indicator of population fitness yet is not typically included in monitoring programs for sensitive species. Additionally, most programs monitor populations at one scale, which can lead to potential mismatches with ecological processes critical to species' conservation. Recently developed methods generating hierarchically nested population units (i.e., clusters of varying scales) for greater sage-grouse (Centrocercus urophasianus) have identified population trend declines across spatiotemporal scales to help managers target areas for conservation. The same clusters used as a proxy for spatial scale can alert managers to local units (i.e., neighborhood-scale) with low genetic diversity, further facilitating identification of management targets. We developed a genetic warning system utilizing previously developed hierarchical population units to identify management-relevant areas with low genetic diversity within the greater sage-grouse range. Within this warning system we characterized conservation concern thresholds based on values of genetic diversity and developed a statistical model for microsatellite data to robustly estimate these values for hierarchically nested populations. We found that 41 of 224 neighborhood-scale clusters had low genetic diversity, 23 of which were coupled with documented local population trend decline. We also found evidence of cross-scale low genetic diversity in the small and isolated Washington population, unlikely to be reversed through typical local management actions alone. The combination of low genetic diversity and a declining population suggests relatively high conservation concern. Our findings could further facilitate conservation action prioritization in combination with population trend assessments and (or) local information, and act as a base-line of genetic diversity for future comparison. Importantly, the approach we used is broadly applicable across taxa.

New strategies for characterizing genetic structure in wide-ranging, continuously distributed species: A Greater Sage-grouse case study.

Characterizing genetic structure across a species' range is relevant for management and conservation as it can be used to define population boundaries and quantify connectivity. Wide-ranging species residing in continuously distributed habitat pose substantial challenges for the characterization of genetic structure as many analytical methods used are less effective when isolation by distance is an underlying biological pattern. Here, we illustrate strategies for overcoming these challenges using a species of significant conservation concern, the Greater Sage-grouse (Centrocercus urophasianus), providing a new method to identify centers of genetic differentiation and combining multiple methods to help inform management and conservation strategies for this and other such species. Our objectives were to (1) describe large-scale patterns of population genetic structure and gene flow and (2) to characterize genetic subpopulation centers across the range of Greater Sage-grouse. Samples from 2,134 individuals were genotyped at 15 microsatellite loci. Using standard STRUCTURE and spatial principal components analyses, we found evidence for four or six areas of large-scale genetic differentiation and, following our novel method, 12 subpopulation centers of differentiation. Gene flow was greater, and differentiation reduced in areas of contiguous habitat (eastern Montana, most of Wyoming, much of Oregon, Nevada, and parts of Idaho). As expected, areas of fragmented habitat such as in Utah (with 6 subpopulation centers) exhibited the greatest genetic differentiation and lowest effective migration. The subpopulation centers defined here could be monitored to maintain genetic diversity and connectivity with other subpopulation centers. Many areas outside subpopulation centers are contact zones where different genetic groups converge and could be priorities for maintaining overall connectivity. Our novel method and process of leveraging multiple different analyses to find common genetic patterns provides a path forward to characterizing genetic structure in wide-ranging, continuously distributed species.

Historical effective population size of North American hoary bat (Lasiurus cinereus) and challenges to estimating trends in contemporary effective breeding population size from archived samples.

BACKGROUND: Hoary bats (Lasiurus cinereus) are among the bat species most commonly killed by wind turbine strikes in the midwestern United States. The impact of this mortality on species census size is not understood, due in part to the difficulty of estimating population size for this highly migratory and elusive species. Genetic effective population size (Ne) could provide an index of changing census population size if other factors affecting Ne are stable. METHODS: We used the NeEstimator package to derive effective breeding population size (Nb) estimates for two temporally spaced cohorts: 93 hoary bats collected in 2009-2010 and an additional 93 collected in 2017-2018. We sequenced restriction-site associated polymorphisms and generated a de novo genome assembly to guide the removal of sex-linked and multi-copy loci, as well as identify physically linked markers. RESULTS: Analysis of the reference genome with psmc suggested at least a doubling of Ne in the last 100,000 years, likely exceeding Ne = 10,000 in the Holocene. Allele and genotype frequency analyses confirmed that the two cohorts were comparable, although some samples had unusually high or low observed heterozygosities. Additionally, the older cohort had lower mean coverage and greater variability in coverage, and batch effects of sampling locality were observed that were consistent with sample degradation. We therefore excluded samples with low coverage or outlier heterozygosity, as well as loci with sequence coverage far from the mode value, from the final data set. Prior to excluding these outliers, contemporary Nb estimates were significantly higher in the more recent cohort, but this finding was driven by high values for the 2018 sample year and low values for all other years. In the reduced data set, Nb did not differ significantly between cohorts. We found base substitutions to be strongly biased toward cytosine to thymine or the complement, and further partitioning loci by substitution type had a strong effect on Nb estimates. Minor allele frequency and base quality bias thresholds also had strong effects on Nb estimates. Instability of Nb with respect to common data filtering parameters and empirically identified factors prevented robust comparison of the two cohorts. Given that confidence intervals frequently included infinity as the stringency of data filtering increased, contemporary trends in Nb of North American hoary bats may not be tractable with the linkage disequilibrium method, at least using the protocol employed here.

Composition and distribution of fish environmental DNA in an Adirondack watershed.

BACKGROUND: Environmental DNA (eDNA) surveys are appealing options for monitoring aquatic biodiversity. While factors affecting eDNA persistence, capture and amplification have been heavily studied, watershed-scale surveys of fish communities and our confidence in such need further exploration. METHODS: We characterized fish eDNA compositions using rapid, low-volume filtering with replicate and control samples scaled for a single Illumina MiSeq flow cell, using the mitochondrial 12S ribosomal RNA locus for taxonomic profiling. Our goals were to determine: (1) spatiotemporal variation in eDNA abundance, (2) the filtrate needed to achieve strong sequencing libraries, (3) the taxonomic resolution of 12S ribosomal sequences in the study environment, (4) the portion of the expected fish community detectable by 12S sequencing, (5) biases in species recovery, (6) correlations between eDNA compositions and catch per unit effort (CPUE) and (7) the extent that eDNA profiles reflect major watershed features. Our bioinformatic approach included (1) estimation of sequencing error from unambiguous mappings and simulation of taxonomic assignment error under various mapping criteria; (2) binning of species based on inferred assignment error rather than by taxonomic rank; and (3) visualization of mismatch distributions to facilitate discovery of distinct haplotypes attributed to the same reference. Our approach was implemented within the St. Regis River, NY, USA, which supports tribal and recreational fisheries and has been a target of restoration activities. We used a large record of St. Regis-specific observations to validate our assignments. RESULTS: We found that 300 mL drawn through 25-mm cellulose nitrate filters yielded greater than 5 ng/µL DNA at most sites in summer, which was an approximate threshold for generating strong sequencing libraries in our hands. Using inferred sequence error rates, we binned 12S references for 110 species on a state checklist into 85 single-species bins and seven multispecies bins. Of 48 bins observed by capture survey in the St. Regis, we detected eDNA consistent with 40, with an additional four detections flagged as potential contaminants. Sixteen unobserved species detected by eDNA ranged from plausible to implausible based on distributional data, whereas six observed species had no 12S reference sequence. Summed log-ratio compositions of eDNA-detected taxa correlated with log(CPUE) (Pearson's R = 0.655, P < 0.001). Shifts in eDNA composition of several taxa and a genotypic shift in channel catfish (Ictalurus punctatus) coincided with the Hogansburg Dam, NY, USA. In summary, a simple filtering apparatus operated by field crews without prior expertise gave useful summaries of eDNA composition with minimal evidence of field contamination. 12S sequencing achieved useful taxonomic resolution despite the short marker length, and data exploration with standard bioinformatic tools clarified taxonomic uncertainty and sources of error.

Development of microsatellite loci for two New World vultures (Cathartidae).

OBJECTIVE: Use next-generation sequencing to develop microsatellite loci that will provide the variability necessary for studies of genetic diversity and population connectivity of two New World vulture species. RESULTS: We characterized 11 microsatellite loci for black vultures (Coragyps atratus) and 14 loci for turkey vultures (Cathartes aura). These microsatellite loci were grouped into 3 multiplex panels for each species. The number of alleles among black vulture samples ranged from 2 to 11, and 3 to 48 among turkey vulture samples.

Long-term persistence of horse fecal DNA in the environment makes equids particularly good candidates for noninvasive sampling.

Fecal DNA collected noninvasively can provide valuable information about genetic and ecological characteristics. This approach has rarely been used for equids, despite the need for conservation of endangered species and management of abundant feral populations. We examined factors affecting the efficacy of using equid fecal samples for conservation genetics. First, we evaluated two fecal collection methods (paper bag vs. ethanol). Then, we investigated how time since deposition and month of collection impacted microsatellite amplification success and genotyping errors. Between May and November 2014, we collected feral horse fecal samples of known age each month in a feral horse Herd Management Area in western Colorado and documented deterioration in the field with photographs. Samples collected and dried in paper bags had significantly higher amplification rates than those collected and stored in ethanol. There was little difference in the number of loci that amplified per sample between fresh fecal piles and those that had been exposed to the environment for up to 2 months (in samples collected in paper bags). After 2 months of exposure, amplification success declined. When comparing fresh (0-2 months) and old (3-6 months) fecal piles, samples from fresh piles had more matching genotypes across samples, better amplification success and less allelic dropout. Samples defecated during the summer and collected within 2 months of deposition had highest number of genotypes matching among samples, and lowest rates of amplification failure and allelic dropout. Due to the digestive system and amount of fecal material produced by equids, as well as their occurrence in arid ecosystems, we suggest that they are particularly good candidates for noninvasive sampling using fecal DNA.

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach.

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse (Centrocercus urophasianus) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Environmental DNA (eDNA) sampling improves occurrence and detection estimates of invasive burmese pythons.

Environmental DNA (eDNA) methods are used to detect DNA that is shed into the aquatic environment by cryptic or low density species. Applied in eDNA studies, occupancy models can be used to estimate occurrence and detection probabilities and thereby account for imperfect detection. However, occupancy terminology has been applied inconsistently in eDNA studies, and many have calculated occurrence probabilities while not considering the effects of imperfect detection. Low detection of invasive giant constrictors using visual surveys and traps has hampered the estimation of occupancy and detection estimates needed for population management in southern Florida, USA. Giant constrictor snakes pose a threat to native species and the ecological restoration of the Florida Everglades. To assist with detection, we developed species-specific eDNA assays using quantitative PCR (qPCR) for the Burmese python (Python molurus bivittatus), Northern African python (P. sebae), boa constrictor (Boa constrictor), and the green (Eunectes murinus) and yellow anaconda (E. notaeus). Burmese pythons, Northern African pythons, and boa constrictors are established and reproducing, while the green and yellow anaconda have the potential to become established. We validated the python and boa constrictor assays using laboratory trials and tested all species in 21 field locations distributed in eight southern Florida regions. Burmese python eDNA was detected in 37 of 63 field sampling events; however, the other species were not detected. Although eDNA was heterogeneously distributed in the environment, occupancy models were able to provide the first estimates of detection probabilities, which were greater than 91%. Burmese python eDNA was detected along the leading northern edge of the known population boundary. The development of informative detection tools and eDNA occupancy models can improve conservation efforts in southern Florida and support more extensive studies of invasive constrictors. Generic sampling design and terminology are proposed to standardize and clarify interpretations of eDNA-based occupancy models.

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies.

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (∼3.5-5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus, and the Clark's Nutcracker, Nucifraga columbiana). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample--a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Landscape genetics of fishers (Martes pennanti) in the Northeast: dispersal barriers and historical influences.

Habitat fragmentation and overtrapping are thought to have resulted in severe population declines for fisher (Martes pennanti) across the northeastern United States, and by the end of the 1930s only 3 remnant populations remained. Subsequent trapping cessation, extensive reintroduction programs, and natural recolonization have helped fishers to reclaim much of their historical range. The degree to which these processes have impacted genetic structure in this species, however, remains unknown. We used 11 microsatellites from tissue samples (n = 432) of fishers to characterize contemporary population structure in light of historical population structure and thus to determine the relative influence of anthropogenic disturbances and natural landscape features in shaping genetic structure of the contemporary population. Our results indicated that 3 well-differentiated contemporary populations are present that correspond well with what would be expected based on their reported history. A course barrier to dispersal appears in the western portion of the study area associated with several lakes including Lake George and Great Sacandaga Lake. Large-scale reintroduction efforts and natural recolonizations have largely had predictable impacts on population structure. An important exception is the substantial impact of the reintroduction of fishers to Vermont.

Species-wide phylogeography of North American mule deer (Odocoileus hemionus): cryptic glacial refugia and postglacial recolonization.

Quaternary climatic oscillations greatly influenced the present-day population genetic structure of animals and plants. For species with high dispersal and reproductive potential, phylogeographic patterns resulting from historical processes can be cryptic, overshadowed by contemporary processes. Here we report a study of the phylogeography of Odocoileus hemionus, a large, vagile ungulate common throughout western North America. We examined sequence variation of mitochondrial DNA (control region and cytochrome b) within and among 70 natural populations across the entire range of the species. Among the 1766 individual animals surveyed, we recovered 496 haplotypes. Although fine-scale phylogenetic structure was weakly resolved using phylogenetic methods, network analysis clearly revealed the presence of 12 distinct haplogroups. The spatial distribution of haplogroups showed a strong genetic discontinuity between the two morphological types of O. hemionus, mule deer and black-tailed deer, east and west of the Cascade Mountains in the Pacific Northwest. Within the mule deer lineage, we identified several haplogroups that expanded before or during the Last Glacial Maximum, suggesting that mule deer persisted in multiple refugia south of the ice sheets. Patterns of genetic diversity within the black-tailed deer lineage suggest a single refugium along the Pacific Northwest coast, and refute the hypothesis that black-tailed deer persisted in one or more northern refugia. Our data suggest that black-tailed deer recolonized areas in accordance with the pattern of glacial retreat, with initial recolonization northward along a coastal route and secondary recolonization inland.

Development and characterization of 12 polymorphic microsatellite loci in the American dog tick (Dermacentor variabilis).

Twelve polymorphic microsatellite loci were developed for the American dog tick (Dermacentor variabilis), an important vector of infectious diseases in humans and animals. Four multiplexed panels comprising the loci were developed and 45 ticks collected from two raccoons (Procyon lotor) were genotyped. The number of alleles per locus ranged from nine to 30, and single locus heterozygosities ranged from 0.18 to 0.93. Data generated using these markers will further our understanding of factors affecting gene flow in D. variabilis, thus helping to elucidate the transmission dynamics of diseases associated with this vector.

Development and characterization of 14 polymorphic microsatellite loci in the raccoon tick (Ixodes texanus).

Fourteen polymorphic microsatellite loci were developed for the raccoon tick (Ixodes texanus). Three multiplexed panels comprising the loci were developed and 50 ticks collected from two infected raccoons (Procyon lotor) were genotyped. The number of alleles per locus ranged from five to 22, and single locus heterozygosities ranged from 0.46 to 0.94. Future research will further our understanding of I. texanus biology and help in elucidating the effects of life-history variation on parasite population genetic structure, using I. texanus as a model organism.

Identification of 24 polymorphic microsatellite markers for the double-crested cormorant (Phalacrocorax auritus).

Twenty-four polymorphic microsatellite markers were developed for the double-crested cormorant (Phalacrocorax auritus). The number of alleles ranged from two to 13 and observed heterozygosities ranged from 0.032 to 0.871. The use of these loci should enable researchers and biologists to learn more about the population structure and ecology of this species.

Isolation of 21 polymorphic microsatellite markers for the Virginia opossum (Didelphis virginiana).

Twenty-one polymorphic microsatellite markers were developed for the Virginia opossum (Didelphis virginiana). The number of alleles ranged from two to 13 and observed heterozygosities ranged from 0.464 to 0.964. Significant heterozygote deficiencies were observed at three loci and null alleles were detected at five loci. Evidence for gametic disequilibrium was observed between three sets of paired loci after a sequential Bonferroni correction was applied. These markers will enable us to investigate the mating tactics, movement behaviour and social structure of Virginia opossum populations inhabiting fragmented agricultural landscapes.